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ABSTRACT 

The radio continuum emission from the Galaxy has a rich mix of thermal and non-thermal emission. 
This very richness makes their interpretation challenging since the low radio opacity means that a 
radio image represents the sum of all emission regions along the line of sight. These challenges make 
the existing narrow-band radio surveys of the Galactic plane difficult to interpret: e.g., a small region 
of emission might be a supernova remnant (SNR) or an Hll region, or a complex combination of both. 
Instantaneous wide bandwidth radio observations in combination with the capability for high resolu- 
tion spectral index mapping, can be directly used to disentangle these effects. Here we demonstrate 
simultaneous continuum and specitral index imaging capability at the full continuum sensitivity and 
resolution using newly developed wide-band wide-field imaging algorithms. Observations were done in 
the L- and C-Band with a total bandwidth of 1 and 2 GHz respectively. We present preliminary results 
in the form of a full-field continuum image covering the wide-band sensitivity pattern of the EVLA 
centered on a large but poorly studied SNR (G55.7 + 3.4) and relatively narrower field continuum and 
spectral index maps of three fields containing SNR and diffused thermal emission. Wc demonstrate 
that spatially resolved spectral index maps differentiate regions with emission of different physical 
origin (spectral index variation across composite SNRs and separation of thermal and non-thermal 
emission), superimposed along the line of sight. The wide-field image centered on the SNR G55.7 -1-3.4 
also demonstrates the excellent wide-field wide-band imaging capability of the EVLA. 
Subject headings: ISM: supernova remnants — techniques: interferometric 



1. INTRODUCTION 

In order to test the wide-band imaging capabilities of 
the EVLA^, we have observed several known and pos- 
sible Galactic supernova remnants (SNRs). There are 
274 Galactic SNRs catalogued (Green 2009), the ma- 
jority of which are classified as 'shell' remnants, show- 
ing limb-brightened synchrotron radio emission, typically 
with S cx Some SNRs are classed as 'filled-centcr' 

remnants (or 'plerions'), as they show emission at radio 
wavelengths that is brightest in the center, due to the 
power input from a central pulsar. These filled-center 
remnants generally have much flatter radio spectra at 
GHz- frequencies, with S oc v' 0.0-0.3^ Finally, there 
are some remnants that are classed as composite rem- 
nants, as they show a mixture of limb-brightened shell- 
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like emission, with a fiUed-center-like core (i.e., a pulsar 
wind nebula) showing a much flatter spectrum. In addi- 
tion to the known remnants, there are many other pos- 
sible and probable Galactic SNRs that have been pro- 
posed, for which additional observations arc needed to 
confirm their nature. In many regions of the Galactic 
plane, particularly near to 6 = 0° and ? = 0° , radio emis- 
sion from SNRs is easily confused with thermal radio 
emission (typicallj^ with S oc v' '^'^) from Hll regions. 
The wide-band of the EVLA allows a variety of spectral 
studies of SNRs to be made, in order to: (1) investigate 
spectral variations across the face of shell remnants, in 
order to study the differences in the shock acceleration 
mechanism at work in different regions of the remnants; 
(2) study the flat spectrum cores of composite remnants; 
and (3) clarify the nature of proposed possible SNRs, 
by disentangling the non-thermal synchrotron emission 
from confusing thermal emission. 

We have observed two smaller known remnants 
G16.7H-0.1 and G21.5-0.9. G16.7+0.1 is a 'compos- 
ite' remnant, 3' across (e.g., Helfand et al. 1989). 
G21.5— 0.9 was classed as a filled-center remnant, only 
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~ 1.5' in extent, until a faint X-ray shell ^ 4' in diam- 
eter was identified by Slane ct al. (2000). Since then it 
has been classed as a composite remnant, although a ra- 
dio counterpart to the X-ray shell has not been identified 
(see Bietenholz et al. 2011). G21.5— 0.9 contains a pul- 
sar, which has a high spin-down luminosity (Gupta et al. 
2005). We also observed a field centered at I = 19.6°, 
b = —0.2°, which contains 'high-probability' SNR can- 
didate, G19.66-0.22, identified by Helfand et al. (2006). 
Finally, we observed one larger known but poorly studied 
shell remnant, G55.7-I-3.4, which is ~ 20 x 25 arcmin^ in 
extent (see Goss et al. (1977)). An old, probably unre- 
lated pulsar lies within the northern edge of this remnant. 

2. OBSERVATIONS AND CALIBRATION 

The observations were made using the EVLA in D- 
array. The L-Band observations used the fuU 1 GHz 
bandwidth covering the frequency range of 1 — 2 GHz. 
The C-Band observations used 2 GHz bandwidth in two 
separate ~ 1 GHz bands. The observations parameters 
for all the fields are listed in Table 1. At the time of 
observations, not all antennas in the array were outfit- 
ted with the new L- and C-Band receivers. Due to this 
and other failures during the commissioning phase of the 
EVLA, typically 20 antennas were used for these obser- 
vations. The L-Band data was recorded in 8 or 16 spec- 
tral windows (SPW) at a frequency resolution of 2 and 

1 MHz respectively with an integration time of 1 s. The 
C-Band data was recorded in 16 SPWs at a resolution of 

2 MHz. High time and frequency resolution was required 
to allow effective removal of radio frequency interference 
(RFI), accurate bandpass calibration and corrections for 
wideband group delay present in the data from EVLA 
WIDAR correlator. To reduce the data volume, after 
initial data flagging the data were averaged to 10 s reso- 
lution in time. The standard flux calibrators 3C286 and 
3C147 were observed for flux calibration and the com- 
pact sources J1822-0938 (for G16. 7+0.1, G19.6-0.2 and 
G21.5-0.9) and J1925+2106 (for G55.7-(-3.4) were ob- 
served at intervals of 30 min for phase and bandpass cali- 
bration. Using known frequency dependent model for the 
fiux calibrators (including spectral index) antenna gains 
for each frequency channel across the observed band- 
width were determined for flux and group delay calibra- 
tion. The antenna gains for each SPW were also deter- 
mined using the phase calibraters. In order to account for 
the spectral index of the calibraters, a model including 
the frequency dependent flux was determined using the 
multi-scale multi-frequency synthesis (MS-MFS) imag- 
ing algorithm which simultaneously makes Stokes-/ and 
spectral index maps. The final complex bandpass shape 
was determined using this model and applied to the data 
for the target fields. This calibration procedure allowed 
calibration of time-varying direction independent gains 
without transferring the effects of the frequency depen- 
dent fiux of the calibraters to the target fields. 

3. WIDE-BAND WIDE-FIELD IMAGING 

RFI affected frequency channels were fiagged through- 
out the observing band. Full SPWs covering the fre- 
quency range 1.516 ^ 1.643 GHz had to be dropped due 
to the presence of strong RFI at i^Band. Ten channels 
at the edges of the SPWs were also fiagged due to the 
roll-off of the digital filters in the signal chain. The rest of 



the frequency channels were used for wide-band imaging. 

Conventional image reconstruction techniques for in- 
tcrferometric imaging fundamentally models the bright- 
ness distribution as a collection of scaleless components 
(amplitude per pixel) leaving errors which limit the imag- 
ing performance. Scale-sensitive methods (Bhatnagar & 
Cornwell 2004; CornwcU 2008) significantly reduce the 
magnitude of such errors by explicitly solving for the 
scale size of emission across the field of view. Imag- 
ing performance using wide-band interferometric obser- 
vations of fields with extended emission is additionally 
limited by the fact that not only does the scale of emis- 
sion vary across the field of view, the spectral proper- 
ties also vary across the field. The dominant source of 
these spectral variations is the change of antenna pri- 
mary beam with frequency and spatial dependence of 
the spectral index of the radio emission. Both effects 
result in variation of the strength of emission as a func- 
tion of frequency in a direction dependent manner. Fur- 
thermore, antenna primary beams are typically rotation- 
ally asymmetric which result in time- varying gains due to 
the rotation of the primary beams with parallactic angle 
for El-Az mount antennas. Scale-sensitive imaging re- 
construction algorithms that can also simultaneously ac- 
count for time and frequency-dependent effects are there- 
fore required for wide-band wide-field imaging with the 
EVLA, particularly in the Galactic plane. For this, two 
new algorithms, namely the MS-MFS (Rau 2010; Rau & 
Cornwell submitted,2011) and the A-Projection (Bhat- 
nagar et al. 2008) algorithm, have been developed. The 
MS-MFS algorithm combines the scale-sensitive decon- 
volution with explicit frequency dependent modeling of 
the emission to map the spatial distribution of the spec- 
tral index in a scale-sensitive manner. The A-Projection 
algorithm on the other hand corrects for the time and 
frequency dependence of the antenna primary beams as 
part of the iterative image reconstruction. 

We used the MS-MFS (Rau 2010; Rau & CornweU sub- 
mitted, 2011) algorithm implemented in the CASA^ post- 
processing package, for wide-band imaging. The spectral 
index information derived using the traditional method 
of using images made at different frequencies make spec- 
tral index images at the lowest resolution and at the 
sensitivity of single channel bandwidth. Spectral index 
images made using the MS-MFS algorithm on the other 
hand does not require smoothing the images to the lowest 
resolution and produces the spectral index images at the 
full sensitivity and resolution offered by the data. This 
must be combined with techniques to correct for various 
time variable PB-effects during image deconvolution to 
achieve noise limited imaging performance. While the 
A-Projection algorithm can be used to correct for these 
PB-effects during image deconvolution, for the images 
presented here we used only the time-averaged PBs (see 
Figure 3) to apply correction post deconvolution. Work 
towards a fully integrated wide-field wide-band imaging 
algorithms is in progress. 

4. RESULTS 

4.1. G16. 7+0.1, G21.5-0.9, G19.6-0.2 

Figure 1 shows the Stokes-/ and Spectral Index im- 
ages of two catalogued SNRs (G16.7-I-0.1 and G21.5-0.9) 

^ http://casa.nrao.edu 
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TABLE 1 
List of Observations Parameters 



Field 


Pointing Center 


Date 


Frequency Range 


Calibrators 




(J2000) 




(GHz) 


Flux, Phase 


G16. 7+0.1 


RA=18:20:57 


2010 Aug 12 


4.49-5.45 


J1331+3030 and J0137+331, 




Dec=-14:19:30 




6.89-7.85=^ 


J1822-0938 


G19. 6-0.2 


RA=18:27:38 


2010 Jul 19 


1.00-2.03 


J1331+3030, J1822-0938 




Dec=-ll:56;32 








G21. 5-0.9 


RA=18:33:32 


2010 Aug 12 


4.49-5.45 


J1331+3030 and J0137+331, 




Dec=-10:34:10 




6.89-7.85=^ 


J1822-0938 


G55. 7+3.4 


RA=19:21:40 


201 Aug 23 


1.00-2.03 


J1331+3030 and J0542+498, 




Dec=+21:45:00 






J1925+2106 



^The two ranges correspond to the two 1-GHz wide frequency 
bands covering total of 2 GHz bandwidth at C-Band. 
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Fig. 1. — Stokes-/ (top row) and Spectral Index images (bottom row). The first two columns show C-Band images of the SNRs G16. 7+0.1 
and G21.5— 0.9 respectively. The resolution in these images is 10" and RMS noise is 11/iJy beam~^ and 30^ Jy beam~^ respectively. The 
last column shows the L-Band image of G19.6— 0.2 with a compact source of thermal emission of positive spectral index surrounded by 
non-thermal emission with a negative spectral index. The resolution in this image is 40" and RMS noise of 0.2m, Jy beam~^. The circle 
indicates the half-power point of the antenna primary beam at the reference frequencies of 6.0 GHz for C-Band images and 1.5 GHz for 
the L-Band image. 



(Green 2009) and one candidate SNR (G19.6-0.2). The 
top row shows the Stokes-/ images and the bottom row 
shows the spectral index images, with the Stokes-/ con- 
tours overlaid on the spectral images. Correction for the 
primary beam response raises the noise away from the 
center. In order to show the detailed brightness distribu- 
tion, we chose to show the Stokes-/ images without cor- 
recting for the primary beam response. The integrated 
flux reported below was determined from PB-corrected 
images. The spectral index images have been corrected 
for the average primary beam effects, since the effect of 
primary beam response in the spectral index images is 
stronger. Image of all the fields covering the full wide- 
band sensitivity show many more compact and diffused 
sources. Due to space limitation, here we show central 
portions of these images. 
The first column has the images at a reference fre- 



quency of 6 GHz, of the composite SNR G16. 7-1-0.1 at a 
resolution of 10". The SNR has a flatter spectrum core 
with an average spectral index of —0.54 ± 0.02 {S oc ly") 
surrounded by a relatively steeper spectrum nebula with 
the spectral index ranging from —0.8 to 1.1. The to- 
tal continuum integrated flux is 1.2 Jy with an RMS 
noise of 11 ^ Jy beam~^ . The diffused sources on east of 
SNR are unrelated but real sources of emission. The im- 
ages in the second column are of the fiUed-centered SNR 
G21.5— 0.9, also at the reference frequency of 6 GHz. The 
RMS noise in this image is '-^ 30/i Jy beam^^ and the 
resolution is 10". The spectral index across this SNR 
is uniform across with a value of +0.12 ± 0.03 and in- 
tegrated flux density of 6 Jy. Also visible in the im- 
age is the compact source called the "northern knot" 
about 2' to the north. The long extended feature ex- 
tending towards the north-west direction is part of the 
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separate catalogued SNR G21.6+0.84. The third col- 
umn contains the L-Band images at a reference frequency 
of 1.5 GHz, of the region centered on the Galactic co- 
ordinates / = 19.6°, b = —0.2°. This field contains a can- 
didate SNR G19.66-0.22 (Helfand et al. 2006). Stokes-/ 
image shown in the top row shows weak diffused emis- 
sion surroimding two compact strong sources, making it 
difficult to determine the nature of the objects based on 
Stokes-/ morphology alone. The spectral index image 
however clearly separates sources of thermal and non- 
thermal emission. The spectral index of the diffused 
emission range from —0.4 to —0.5. The spectral index of 
the compact source (G19.67— 0.15) towards the norther 
edge of the diffused nebida has a spectral index of —0.3. 
The superimposed compact thermal source; near the cen- 
ter is the known ultra-compact Hii region G19.61— 0.23 
(e.g.. Wood & Churchwell (1989); Hofner & Churchwell 
(1996)) also with a spectral index of -1-0.3. Based on mor- 
phology of the spectral distribution we conclude that this 
field contains an SNR, possibly of filled-centered type. 

4.2. G55.7+34 

The SNR G55. 7-1-3.4 is classified as an "incomplete- 
shell SNR" with a pulsar within the boundary of the 
remnant based on the only existing radio image at 610 
and 1.415 GHz by Goss et al. (1977). These are the 
only image published in the literature of this SNR made 
over three decades ago using the Westerbork Synthesis 
Radio Telescope. The full 2° x 2° Stokes-/ L-Band im- 
age centered on this SNR at the reference frequency of 
1.5 GHz covering the full wide-band sensitivity pattern 
of the EVLA antenna is shown in Figure 2. This is the 
highest resolution and sensitivity image of this region to 
date. The wide-band sensitivity pattern of the EVLA 
extended significantly beyond even the first sidelobe of 
the PB at the lowest frequency in the band. The time- 
averaged wide-band sensitivity pattern of the EVLA for 
the equivalent frequency range and field-of-view is shown 
in the left panel of Figure 3. The one-dimensional slice 
through this pattern in the right panel, shows the sensi- 
tivity pattern extending significantly beyond the "main- 
lobe" at the level of 5 — 15% of the peak value. With 
sources, some strong, detected up to 1° away from the 
center of the image, imaging the full field was required 
to reach noise-limited imaging at the center of the field. 
Due to this wide field of view, even for D-array observa- 
tions, correction for the effects of non-coplanar baselines 
had to be applied. This was done using the W-Projection 
algorithm (Cornwell et al. 2008). The Stokes-/ image of 
the SNR shows a nearly complete shell-type SNR. At a 
resolution of 30", this image for the first time, shows a 
network of filamentary structures along with smoother 
diffused emission filling the boundary of the shell. The 
spectral index image of this SNR however was not re- 
liable at all scales since the large scale emission is not 
well constrained by the shortest measured baselines at 
all frequencies. The filaments however show up as re- 
gions of steeper spectral index compared to an aver- 
age spectral index of —0.6 of the surrounding regions. 
The strong compact source (G55.78+3.5) on the north- 
ern edge of the shell in the Stokes-/ image is an older 
unrelated pulsar (PSR J1921+2153). In the spectral in- 
dex map (not shown here), this pulsar shows up as a 
steep spectrum compact source with a spectral index of 



—2.2. The integrated flux of the shell, including the pul- 
sar and other foreground or background sources super- 
imposed across the shell is 1.0 Jy. The rms noise in the 
image is 10 fiJy beam~^. 

5. DISCUSSION 

The technical goals of this pilot project were to assess 
our low-frequency wide-band imaging capabilities for tar- 
gets in the Galactic plane. At i-Band and C-band, the 
Galaxy is replete with compact and diffuse emission of 
both thermal and non-thermal origin. Sources often have 
a mixture of thermal and non-thermal emission, resulting 
in spectral indices that vary with position across a given 
source and with spatial scale. Conventional imaging al- 
gorithms that rely on multiple narrow-band images to de- 
rive spectral-index information are limited to the angular 
resolution allowed by the lowest measured frequency, and 
the single-channel sensitivity. The MS-MFS algorithm 
was designed to optimally use the multi-frequency mea- 
surements to do spectral-index mapping across diffused 
sources at the highest sensitivity and angular resolution 
offered by the data. 

The following is a summary of our current status, based 
on the imaging results shown in the previous section. 
Our targets were chosen based on existing information 
about their spatial and spectral structure, in order to 
test the analysis methods in preparation for a Galactic 
plane survey. 

1. Separation of thermal versus non-thermal; steep 
versus steeper spectra : we see clear morphologi- 
cal separations between regions expected to have 
different spectral structure from the surroundings. 
These differences range from Aa « 0.6 between 
thermal and non-thermal regions, down to Aa « 
0.2 between structures thought to be cores and 
shells. In regions with high S/N ratios (> 100), 
the estimated uncertainty in the spectral index is 
±0.03. In regions with low signal-to-noise (< 100), 
the uncertainty in the spectral index rises rapidly, 
and all the images presented here used a Stokes-/ 
S/N threshold of 100. 

2. Image-reconstruction errors : errors in the spectral 
reconstructions are dominated by multi-scale de- 
convolution artifacts that decrease with appropri- 
ate choices of multi-scale imaging parameters (a set 
of scale-sizes with which to model the spatial struc- 
ture). Error estimates from Monte-Carlo simula- 
tions (by varying the choice of spatial scales) range 
from ±0.02 when scales are chosen appropriately, 
up to 0.2 in the extreme case where diffuse emission 
is modeled by a series of ^-functions. These sim- 
ulations were done for the EVLA D-configuration 
and the errors depend on the precise uv-coverage 
just as they depend for the standard imaging algo- 
rithms. The MS-MFS algorithm reconstructs the 
frequency dependence of the flux by a polynomial 
fit across the frequency range on a per component 
basis. The errors on the estimated spectral index 
therefore also depend on the errors of the fitted 
coeffcients. These errors depend on variations in 
S/N ratio (across the frequency range or due to 
the scale of the emission) in the same way as the 
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Fig. 2. — Stokes-/ image at a reference frequency of 1.5 GHz. Data in the frequency range 1.256-1.905 GHz were included for this image. 
For this frequency range, the HPBW of the primary beam changes from 32' to 21'. This 2° X 2° image covers the full-field corresponding 
to the wide-band sensitivity pattern of the EVLA (see Figure 3) of the field containing the SNR G55.7-|-3.4. 



errors on the coefficients of any least-squares poly- 
nomial fitting algorithm (sec Press ct al. (2007); 
Rau & CornwcU (submitted, 2011) for further de- 
tails). Further, at the largest spatial scales sampled 
by the interferometer (which, of course, is telescope 
dependent), the spectrum is unconstrained by the 
data and the MS-MFS model does not provide ad- 
equate additional constraints. Out of our exam- 
ples, G55.7+3.4 and G16. 7+0.1 contain structme 
at spatial scales that fall within the unsamplcd cen- 
tral region of the interferometer sampling function, 
aleading to overall systematic errors in the spectral 
index across the source. These errors are limited to 
the largest scales, but in the current implementa- 
tion of the algorithm, it is not possible to analyse 
the scales separately. In our examples, G55. 7+3.4 



is affected significantly by this error and spectra of 

only the most compact emission can be trusted. Fi- 
nally, in addition to errors in the spectral index, off- 
source artifacts arc present in the Stokes-/ images 
beyond dynamic-ranges of a few thousand (mea- 
sured as the ratio of the peak brightness to the 
peak residual near the source). Dynamic-ranges 
for the images presented here range from 2000 for 
G16. 7+0.1 to 10000 for G19. 6-0.2. In most cases, 
wide-band self-calibration was required to achieve 
these limits. 



3. Wide-field wide-band sensitivity of the instrument, 
and related errors : across a 2:1 bandwidth, the 
size of the antenna power-pattern changes enough 
to produce an average beam that is significantly 
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Fig. 3. — The EVLA time-averaged wide-band sensitivity pattern covering the same frequency range and the area shown in Figure 2. 
The left panel shows the pattern in gray scale with overlaid contours at 50%, 25%, 15%, 10% and 6% of the peak value. The right panel 
shows a one-dimentional cut through the pattern shown on the left. 



different from the beams at individual frequencies. 
The most-prominent difference is the near lack of 
a first null, with sensitivity at the few percent 
level out to almost three times the HPBW. The 
G55.7-I-3.4 field is an example of this wide- field sen- 
sitivity pattern (the reference-frequency primary 
beam has a HPBW of 30'). For sources far from 
the phase-center, the w-term causes the dominant 
wide-field error, and these images (in particular 
G55.7-I-3.4) are a result of W-Projection combined 
with MS-MFS. The next dominant error (in Stokes- 
I and spectral index) is due to the frequency depen- 
dence of the primary beam. All the spectral-index 
maps shown here have been corrected for this ef- 
fect using the time-averaged PB-model, and can 
be trusted out to the 25% point of the reference- 
frequency primary-beam. The time-variability of 
the primary-beams has not been accounted for, and 
is thought to be the cause of the next dominant er- 
ror (visible as artifacts around bright sources about 



1° from the phase center in G55.7+3.4). 

The next steps in our analysis include combining multi- 
configuration and multi-band data for the G16. 7-1-0.1 
and G19.6— 0.2 fields to introduce more constraints dur- 
ing image-reconstruction, to complete the software inte- 
gration required to combine MS-MFS and W-Projection 
with the A-Projection algorithm to correct for time- 
varying primary-beams, and further tests with mosaick- 
ing observations and how the additional information aids 
the image reconstruction. These steps are in accordance 
with our long-term goal of doing a wide-band mosaiced 
survey of the Galactic plane at low frequencies. 



We used the ATNF Pulsar Catalogue (Manch- 
ester et al. 2005), accessible via the URL 
http : //www. atnf . csiro . au/people/pulsar/psrcat/ 
for the co-ordinates of the pulsar in the G55.7-I-3.4 field. 
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